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ABSTRACT 

A new, complete, theoretical rotational and vibrational line list for the 
A 2 n <— X 2 X + electronic transition in MgH is presented. The list includes 
transition energies and oscillator strengths for all possible allowed transitions and 
was computed using the best available theoretical potential energies and dipole 
transition moment function with the former adjusted to account for experimental 
data. The A <— X line list, as well as new line lists for the B' 2 £ + <— X 2 £ + 
and the X 2 E + <— X 2 S + (pure rovibrational) transitions, were included in 
comprehensive stellar atmosphere models for M, L, and T dwarfs and solar-type 
stars. The resulting spectra, when compared to models lacking MgH, show that 
MgH provides significant opacity in the visible between 4400 and 5600A . Fur- 
ther, comparison of the spectra obtained with the current line list to spectra 
obtained using the line list constructed by Kurucz (1993a) show that the Kurucz 
list significantly overestimates the opacity due to MgH particularly for the bands 
near 5150 and 4800A with the discrepancy increasing with decreasing effective 
temperature. 
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1. Introduction 

The study of the spectra of cool stars requires detailed knowledge of molecular opacities. 
This includes important absorbers such as TiO, CO, and water vapor, which have bands that 
cover large wavelength ranges and are very important for the structure of the atmosphere 
due to their overall cooling or heating effects. 

In addition, there are a number of molecules that have bands covering comparatively 
small wavelength ranges (e.g., a few 10 or 100 A). Many of them are trace molecules that 
have only small effects on the overall physical conditions inside the atmosphere but that are 
important for spectral classification and for the determination of stellar parameters such as 
effective temperatures, gravities and abundances. Unfortunately, important molecular data 
such as energy levels, bound-free, and bound-bound cross-sections are only poorly known or 
not known at all for a number of these trace molecules. We have therefore started a project 
to update or provide for the first time molecular data of astrophysical interest for important 
trace molecules and consider in this work MgH. These data will be computed using state-of- 
the-art molecular physics codes and should improve our ability to model and analyze cool 
stellar atmospheres considerably. 

It is important to assess the quality of the computed molecular data. This is best done 
by comparing to experimental results; however, this is only possible for very few molecules 
of astrophysical interest. In addition, in many cases the temperature range of astrophysical 
importance is higher than what can be reached with current experimental setups. Therefore, 
indirect methods of testing and evaluating the molecular data are useful. In this paper, we 
use the general-purpose stellar atmosphere code PHOENIX to calculate model atmospheres 
and synthetic spectra with and without the new molecular data. The results of these calcu- 
lations can then be used to assess the importance of particular molecular opacities on the 
structure of the atmosphere. The synthetic spectra can be used to verify the correct strength 
of the computed bands when compared to observational data. This procedure introduces 
uncertainties such as the treatment of the equation of state (e.g., the molecular data used in 
it), the treatment of lines and line profiles and the assumed parameters of the comparison 
star. However, differential analyses circumvent many of these problems and should allow a 
reasonable evaluation of the molecular bound-free and free-free data. 

The electronic bands of magnesium hydride have been detected over a wide range of 
stellar atmospheres including the photosphere of the sun (Sotirovski 1972), sunspot umbrae 
(Wallace et al. 1999), F-K giants in the Milky Way halo and halos of other Local Group 
galaxies (Majewski et al. 2000), and nearby L-dwarfs (Reid et al. 2000). MgH lines can be 
used as indicators of surface gravity in late-type stars (Bonnell & Bell 1993) and to determine 
magnesium isotope abundances (Wallace et al. 1999; Gay & Lambert 2000). 
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The spectrum of MgH has been extensively studied in the laboratory for many decades 
(Balfour & Cartwright 1976; Balfour & Lindgren 1978; Bernath et al. 1985; Wallace et al. 
1999, and references therein) and has received some theoretical attention (Saxon et al. 1978; 
Kirby et al. 1979, and references therein). However, modern stellar atmosphere calculations 
require extensive, and complete, molecular line lists as molecular band absorption is the pri- 
mary source of line-blanketing in cool stellar atmospheres, particularly M dwarfs. For many 
molecules, including MgH, the only source for complete line lists is the extensive compila- 
tions of Kurucz (1993a). While these compilations are highly valuable to stellar modelers, 
the methods necessary to compute 100s of millions of lines require a number of approxima- 
tions which at times are severe. In this work, we apply fully quantum-mechanical techniques 
to compute the complete line list for the A 2 U <— X 2 £ + transition of MgH. The parameters 
of the calculation are adjusted to force agreement with available experiments. However, our 
goal is to reproduce the global MgH opacity, hence we cannot claim spectroscopic accuracy 
for a particular line. The A <— X line list constructed in this work is combined with line 
lists for the B' 2 £ + <— X 2 E + and the X 2 £ + <— X 2 S + transitions computed by Skory et 
al. (2002) and tested in a range of stellar atmosphere models. An overview of the theory of 
molecular rotational lines is presented in section 2 with the results of the line list calculations 
and stellar models are given in section 3. We present our conclusions in section 4. 



2. Molecular theory 

2.1. Potential curves and dipole transition moments 

For the purpose of the present calculations, accurate ab initio potential-energy surfaces 
given by Saxon et al. (1978) have been used for both the A 2 U and X 2 S + electronic 
states of MgH. Their configuration interaction treatment includes singly- to triply-excited 
configurations (SDTCI) in a large Slater-type basis set, leading to good agreement of the 
derived spectroscopic constants with experimental data. 

In order to bring the calculated potential curves into better agreement with experiment, 
shifts of +8.740 x 10~ 4 a.u. 1 and —3.3665 x 10~ 3 a.u. have been applied in the present work to 
the Saxon et al. (1978) X 2 S + and A 2 U energies, respectively. Details about this procedure 
are given in a separate publication for the B' <— X electronic transition in MgH (Skory et 
al. 2002). In this way, the dissociation energies Dq coincide with the experimental values 
given by Balfour & Lindgren (1978) for the ground state and by Balfour & Cartwright (1976) 



Atomic units are used throughout this section unless otherwise stated. 
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for the A 2 II state. The dissociation energies adopted here are 10,243.26 and and 12,903.71 
cm -1 for the X 2 E + and A 2 U of 24 MgH, respectively. The relative energies between the two 
potential curves were further shifted to match the energy difference corresponding to the pure 
(0, 0) vibronic transition. The value of T = 19278.13 cm -1 , which corresponds to the average 
of the Qi(0) and Q2(l) transition energies measured by Balfour & Cartwright (1976), was 
adopted. We note that this value for T does not agree with the more recent measurements 
of Bernath, Black, & Brault (1985) who find T = 19284.65 cur 1 , though it is unclear how 
their value was arrived at. In fact, their assumption of a Hund's coupling case (a) is in 
obvious contradiction with the small and regular A-doubling parameters, po = 0.0258 cm -1 
and q = 0.00178 cm -1 . Such a case for a 2 n <— 2 S + transition is more typically characterized 
by Hund's case (b) (Herzberg 1950, Figures 122 and 124). Further, (Bernath et al. 1985) 
suggest that a number of low-J lines, presumably -Pi(l), -^2(1), ^2(0), Q2(l), ^2(2), and 
P 2 (2), measured by Balfour & Cartwright (1976) do not exist apparently assuming Hund's 
case (a). However, we find that under Hund's case (b) all but -Pi(l) and Q 2 (0) do in fact 
exist. We therefore adopted the value of T deduced from Balfour & Cartwright (1976). 

Out of the range of internuclear separations R = 2.2 a G to 9.5 a considered by Saxon 
et al. (1978), the potential curves have been extrapolated in two different ways. On the one 
hand, for internuclear distances R > 9.5 a G , a smooth fit to the ab initio potentials has been 
performed using the average long-range interaction potential 

where Cq, C§ and C10 are the usual van der Waals coefficients corresponding to the dipole- 
dipole, dipole-quadrupole, and the sum of dipole-quadrupole and dipole-octupole interac- 
tions, respectively. For the X 2 S + ground state, the coefficients described in Skory et al. 
(2002) were adopted as summarized in Table 1. Since to our knowledge, no data have been 
reported for the van der Waals constants of the A 2 Ii state of MgH, they were estimated 
using a technique based on the London formula as described for the B' 2 S + state by Skory et 
al. (2002). On the other hand, for internuclear distances R < 2.2 a G , the potential curves of 
both the X 2 S + and the A 2 U electronic states have been fitted to the short-range interaction 
exponential form Aexp(-BR) + C . 

In a similar way, the dipole transition moment for A 2 n <— X 2 S + and the dipole moment 
of the X state calculated by Saxon et al. (1978) have been used over the range R = 2.2 a G to 
9.5 a D , and extrapolated by an exponential fit for both the short- and long-range interactions. 
The dipole transition moment and the dipole moment were both smoothly forced to zero at 
the united-atom and separated-atom limits. 
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2.2. Rotational and band oscillator strengths 

Throughout the present study, we have adopted the now well established point of view 
expressed by Whiting & Nicholls (1974) with respect to the preferred way in which dipole 
matrix elements of a rotational transition should be separated into rotational and electronic- 
vibrational parts. 

For the A 2 IT <— X 2 S + absorption band system of the present study, the rotational 
oscillator strength can be expressed as given by Larsson (1983), 

fab _ 2 Ap Sj'{J") | n AX |2 / \ 

Jv'J',v"J" - -AiVj'yj"^,, - - \-L>v'J',v"J"\ > \ Z ) 

where D AX , v ,,j„ =< Xv> J'\D AX (R)\xv" J" > is the rovibrational matrix element of the electric 
dipole transition moment D AX (R) responsible for absorption from the X 2 S + into the A 2 U 
electronic state and Xvj are the rovibrational wave functions. The Honl-London factors 
Sj/(J") are defined according to Whiting & Nicholls (1974) as 

r (J" - l)/2, J' = J" - 1 (P-branch) 
Sj,{J") = I (2J"+l)/2, J' = J" (Q-branch) (3) 
1 (J" + 2)/2, J' = J" + 1 (R-branch). 



For the sake of comparison with the band oscillator strength values calculated by Kirby 
et al. (1979), we used the following relation between band and rotational oscillator strengths, 
as defined in Eq.(2), 



/•'/'/ fab 

Jv'v" - ~ab c I T>'\ Jv'J',v"J" 
9j',J"^J'{ J ) 

(2-<Jo )A « +A 0(2J"+l) 



(2-5o,a»)Sj,(J») 



fv'J',v"J"i (4) 



where gji jn is a degeneracy factor arising from spin-splitting and A-doubling in both final 
and initial electronic states. 



3. Results and discussion 

3.1. The A — X rovibrational line list 

The energy levels together with the rovibrational wavefunctions Xv' J'(R) an d Xv"j"(R) 
of the final and initial electronic state have been calculated by solving the radial nuclear 
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equation by standard Numerov techniques (Cooley 1961). These calculations have been 
performed using a step of 1 x 10~ 3 a.u. for the integration over internuclear distances from 
R = 0.5 a.u. to 200 a.u. For 24 MgH, the reduced mass 0.9671852 u 2 = 1763.064 a.u. (Huber 
& Herzberg 1979) was adopted. 

In Table 2, the vibrational binding energies of the A 2 I1 state calculated for the present 
study are given together with the corresponding values of AG(v' + 1/2) = G(v' + 1) — G(v'). 
The latter are found to be in excellent agreement with the previous calculations of Saxon 
et al. (1978) up to v' = 4. For higher vibrational levels, a maximum difference of 125 cm -1 
occurs for v' = 8. Moreover, our fitting procedure to the long-range interaction potential 
yields an additional vibrational level near threshold to give v' max = 13. The maximum J' for 
each v' is also shown in Table 2 giving a total of 435 rovibrational states. 3 Comparison to 
the limited experimental data shows the calculated AG's to be smaller by up to 45 cm -1 
which suggests that the Saxon et al. (1978) potential is somewhat broader than reality. 

The band oscillator strength values for the band system A 2 II <— X 2 S + are given in 
Table 3 for the transitions between the vibrational states v " = — 3 and v' — — 6. The 
current calculations satisfactorily reproduce the theoretical results of Kirby et al. (1979). 

Figure 1 shows the line absorption rotational oscillator strengths 4 for the A 2 U <— 
X 2 S + transition as a function of the energy of the absorbed photon. It is worth noting 
the presence of a series of peaks which correspond to the J' = 1 <— J" = lines of the 
R-branch [R(0)\ whose intensities are larger by a factor two than all other lines in the same 
wavelength region. A detailed examination of the wavelength distribution of these peaks 
reveals that they are to be classified into three categories, according to whether they are 
Av — v' — v" — +1, or — 1 transitions as shown in Figure 1. Apart from these peaks, 
the R-branch contribution dominates over the wavelength ranges 19, 800 — 20, 300 cm -1 and 
21,050 — 21,300 cm -1 . Other lines observed in Figure 1 are to be assigned mainly to the 
Q-branch transitions, which are about two times more intense than the P-branch over the 
whole wavelength range considered. However, the P-branch appears to play a major role for 
the energy range 19, 170 - 19, 280 cm" 1 and 20, 450 - 20, 690 cm" 1 . 

The product gj>,j"-f v 'j',v"j", the so-called gf- value, is plotted in Figure 2 as a function 
2 In atomic weight units, Aston's scale 

3 A discussion of the number and accuracy of the X 2 S + rovibrational states is given in Skory et al. 
(2002). 

4 Thc complete list of MgH oscillator strength data is available online at the UGA Molecular Opacity 
Project database website http://www.physast.uga.edu/ugamop/ 
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of the absorbed photon energy were comparison is made to the line list of Kurucz (1993a). 
To the best of our knowledge, the Kurucz (1993a) calculations were performed with a model 
rotational Hamiltonian using spectroscopic constants (Kurucz 1993b). Thus, the differences 
arising between both sets of results can be mainly explained by the fact that the model 
Hamiltonian method, though satisfactory for low-lying rovibrational levels, may lose accuracy 
as v and/ or J increase. Further, it seems that the same J max was used for all vibrational levels 
(Kurucz, private communication, 2002) which is clearly not the case as shown in Table 2. 
This appears to explain why the band-series in the Kurucz (1993a) calculations extend to 
larger photon energies. 



3.2. Atmosphere models 

The models used for this work were calculated as described in Allard et al. (2001). These 
models and their comparisons to earlier versions are the subject of a separate publication 
(Allard et al. 2001) and we thus do not repeat the detailed description of the models here. 
However, we will briefly summarize the major physical properties. The models are based 
on the Ames H 2 and TiO line lists by Partridge k Schwenke (1997) and Schwenke (1998) 
and also include as a new addition the line lists for FeH by Phillips & Davis (1993) and for 
VO and CrH by R. Freedman (NASA-Ames, private communication). The models account 
for equilibrium formation of dust and condensates and include grain opacities for 40 species. 
In the following, the models will be referred to as "AMES-dusty" for models in which the 
dust particles stay in the layers in which they have formed and "AMES-cond" for models 
in which the dust particles have sunk below the atmosphere from the layers in which they 
originally formed. We stress that large uncertainties persist in the water opacities for parts 
of the temperature range of this work (Allard et al. 2000). However, almost all MgH bands 
are in the optical and are thus affected minimally by the quality of the water opacities. 

In addition to the opacity sources listed above and in Allard et al. (2001, and references 
therein) we added the new A- X bound-bound radiative transition data from this work and 
the B'-X and X-X data from Skory et al. (2002) for all isotopes 24 MgH, 25 MgH, and 26 MgH 
to our opacity database. In order to assess the effects of the new MgH data, we compare 
spectra calculated with these opacity sources to spectra calculated with the MgH line data 
provided by Kurucz in his list of spectral lines of diatomic molecules (Kurucz 1993a). In 
addition, we have computed a number of models with and without MgH data. The original 
AMES grid was calculated for effective temperatures of M, L, and T dwarfs. The hotter 
models in this work are based on the same physics as the AMES grid and merely differ 
in the effective temperature. The models used in the following discussion were all iterated 
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to convergence for the parameters indicated. The high resolution spectra which have the 
individual opacity sources selected are calculated on top of the models. The MgH bands are 
too localized in a region with little flux, or too weak to influence the temperature structure 
of the atmosphere. 

In Figure 3 we show comparisons between spectra from models using no MgH data, 
spectra from models using the new MgH line list and spectra from models using the Kurucz 
(1993a) MgH line list in the spectral region where the MgH bands are most prominent. The 
comparisons have been done at effective temperatures of 2000 K, 3000 K, and 4000 K to 
sample the temperature range in which MgH is visible in the spectrum. For convenience, we 
only did the comparison for models with log(g)=5.0. As can be seen, the Kurucz (1993a) 
line lists overestimates the opacity due to the inclusion of non-existent levels with high J- 
values and the use of a model Hamiltonian approach. The differences between the spectra 
with the new and Kurucz (1993a) MgH line lists are very similar for the AMES-Cond and 
AMES-Dusty models although the overall flux level and the overall flux shape are different 
for AMES-Cond and AMES-Dusty in the optical (Note that reversing the 1.0 dex offset 
applied in the plot is not sufficient to match the two 2000 K models). Another view is 
presented in figure 4 which shows the relative difference in the spectra. As can be seen, 
there are significant differences spreading among both the A — X and the B' — X transitions. 

Since the MgH data presented by (Kurucz 1993a) was intended for solar-type stars we 
also compared the resulting spectra for an effective temperature of 5800 K and log(g)=4.5, 
relevant parameters for solar-type stars. As can be seen as part of Figure 4, the differences 
are negligible. Although the effect of the non-existent levels with high J-values should be 
more important for hotter temperatures, the decrease in the abundance of MgH at higher 
temperatures makes the effect unobservable. Finally, we want to note that for the 2000 K 
AMES-Cond model the X-X band adds a small amount of opacity at around 4/iin and 10/im. 

4. Conclusion 

Using a combination of theoretical and experimental data on the potential energies 
and dipole transition moment of MgH, a comprehensive theoretical vibrational-rotational 
line list for the A 2 n <— X 2 S + transition was constructed. When using the new A <— X 
line data and the new B' 2 S+ <- X 2 £+ and the X 2 S+ <- X 2 S+ line data of Skory et 
al. (2002) in synthetic spectrum calculations, we find significant differences in the opacity 
when comparing the spectra to calculations using the existing data of Kurucz (1993a). The 
differences are largest for effective temperatures pertaining to L and M type stars and can 
easily be seen in low resolution work. For hotter stars, of K and G type, the differences are 
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less pronounced and high resolution spectra are required to notice the improvements for the 
hottest stars. 
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Fig. 1. — Rotational oscillator strengths f v 'j' jV »j» as a function of wavelength of the absorbed 
photon energy. 
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Fig. 2. — gj',j"-f v 'j',v"j" values as a function of the absorbed photon energy, (a) calculations 
of Kurucz (1993a); (b) present calculations. 
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Fig. 3.— Spectra of the strongest MgH bands for models with T eff =4000 K, T cfr =3000 K, 
T cfr =2000 K (AMES-Cond model) and T cff =2000 K (AMES-Dusty model) (from top to 
bottom). The 2000 K AMES-Cond models have an artificial offset of +1.0 dex for better 
legibility. All models have log(g)=5.0 and are AMES-Cond unless noted otherwise. The black 
spectra are from models without MgH data, the dark grey spectra are from models with the 
new MgH data and the light grey spectra are from models using the Kurucz (1993a) MgH 
data. The resolution of the calculated spectra are much higher. For illustrative purposes the 
spectra have been reduced to a resolution of R=80 at 5000 A. 
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Fig. 4. — Relative flux differences between models when changing from the Kurucz (1993a) 
MgH data to the new MgH data. The models are from top left to bottom right 5800 K, 
4000 K, 3000 K and 2000 K. All models are AMES-Cond and have log(g)=5.0, except the 
5800 K model has log(g)=4.5. The resolution of all spectra have been reduced to R=500 
before calculating the ratio. 
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Table 1. Asymptotic Separated- Atom and United- Atom Limits 





Separated-Atom 










Molecular 


Energy (eV) 








United 


State 


Atomic States Theory Expt 


c 6 


c 8 


Cio 


Atom (Al) 


X 2 £+ 
B' 2 £+ 


Mg(3s 2 1 ,S)+H(ls 2 S) 0.0 0.0 
Mg(3s3p 3 P°)+H(ls 2 S) 2.596 2.714 c 
Mg(3s3p 3 P°)+H(ls 2 S) 2.593 2.714 c 


57.8 a 
56.9 d 
56.9 d 


2490 b 
2451 d 
2451 d 


115000 b 
113205 d 
113205 d 


3s 2 3p 2 P° 
3s 2 3p 2 P° 
3s 2 4s 2 S 



a A. Derevianko (private communication, 2001). 
b From Standard & Certain (1985). 

c Weighted average of the 3 P , 3 Pi, and 3 P 2 term energies from NIST Atomic Spectra 
Database (1999). 

d Estimate, Skory et al. (2002). 
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Table 2. Vibrational binding energies' 1 and AG(v' + 1/2) in cm 1 for the A 2 IT state 



AG(v' + 1/2) 



v' 


V 

"max 


B. E. b 


Theory b 


Saxon c 


Expt. d 





49 


12903.7 


1494.3 


1493.9 


1533.9 


1 


47 


11409.4 


1423.2 


1423.2 


1466.1 


2 


44 


9986.2 


1349.0 


1349.5 


1394.4 


3 


42 


8637.2 


1273.2 


1273.0 




4 


39 


7364.0 


1196.4 


1191.8 




5 


36 


6167.5 


1132.0 


1103.9 




6 


33 


5035.5 


1082.9 


1007.4 




7 


30 


3952.6 


1013.1 


900.2 




8 


27 


2939.5 


905.7 


780.2 




9 


23 


2033.8 


744.7 


645.7 




10 


19 


1289.1 


573.3 


495.0 




11 


15 


715.8 


414.1 


330.9 




12 


11 


301.8 


170.1 






13 


6 


131.7 









a Binding energies are given for rotationless vi- 
brational levels 
b This work 

c Saxon et al. (1978) 
d Balfour & Cartwright (1976) 
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Table 3. Band oscillator strengths* and transitions energies^ for the A 2 U <— X 2 S + Band 

System 





( V2sp + 


^ „,// _ n 
) v — u 




V 












V 


= 3 




„/ ( A2yi\ 
V yA 11 J 


TP 

WO 


f 

Jv'O 




TP 


Jv'l 






f 

Jv'l 




TP 

&v'3 


fv'3 





19290 


1.616( 


-1) 


17867 


8.30{ 


-3) 


16506 


5.2( 


-4) 


15209 


4.1( 


-5) 




19292 


1.61( 


-1) 


17860 


8.61( 


-3) 


16492 


5.3( 


-4) 


15190 


4.1( 


-5) 


1 


20784 


4.88(- 


-3) 


19361 


1.424 


;-i) 


18000 


1.49( 


-2) 


16702 


1.48( 


-3) 




20826 


5.14(- 


-3) 


19394 


1.41( 


-i) 


18025 


1.59( 


-2) 


16724 


1.52( 


-3) 


2 


22207 


1.0( 


-5) 


20784 


8.62{ 


-3) 


19422 


1.246( 


-1) 


18125 


1.99( 


-2) 




22288 


3.9( 


-6) 


20856 


9.38( 


-3) 


19490 


1.24( 


-1) 


18185 


1.99( 


-2) 


3 


23555 


6.9( 


-7) 


22132 


5.9( 


-5) 


20770 


1.11( 


-2) 


19473 


1.082( 


-1) 




23682 


4G 


■7) 


22250 


4.4{ 


-5) 


20884 


1.14( 


-2) 


19579 


1.09( 


-1) 


4 


24828 


1.4(- 


-8) 


23405 


2.2( 


-6) 


22043 


2.0( 


-4) 


20747 


1.25( 


-2) 




25002 


< 1.0(- 


■7) 


23570 


7.0{ 


-6) 


22204 


2.4( 


-4) 


20899 


1.27( 


-2) 


5 


26024 


1.3( 


-9) 


24601 


1.3{ 


-7) 


23240 


2.8( 


-6) 


21943 


5.0( 


-4) 




26243 


< 1.0(- 


■7) 


24811 


< 1.0(- 


-7) 


23445 


5.0( 


-6) 


22140 


5.2( 


-4) 


6 


27156 


1.2( 


-9) 


25733 


8.7{ 


-9) 


24371 


6.5( 


-7) 


23074 


1.2( 


-6) 




27398 


< 1.0(- 


7) 


25966 


4.0( 


-7) 


24600 


2.0( 


-6) 


23295 


1.0( 


-6) 



*Band oscillator strengths are given in a.u. and are calculated for rotational quantum 
numbers J" = and J' — 1. Our results are listed on the first line, with the previous 
calculations of Kirby et al. (1979) below. Notation: x(—n) = x x 10~™ 

''"Energies E v i v ii of the absorbed photon are given in cm" 1 . Our values are listed on the 
first line, with the experimental transition energies given by Balfour & Cartwright (1976) 
below (generally as listed in Kirby et al. (1979)), for J" = and J' = I. 



